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Abstract. I study the dynamical effects due to the Brans-Dicke scalar field at 
the early stages of a supposedly anisotropic Universe expansion in the scalar-tensor 
cosmology of Jordan-Brans-Dicke. This is done considering the behaviour of the 
general solutions for the homogeneous model of Bianchi type VII in the vacuum 
case. I conclude that the Bianchi- VIIo model shows an isotropic expansion and that 
its only physical solution is equivalent to a Friedman-Robertson- Walker spacetime 
whose evolution can, depending on the value of the JBD coupling constant, begin in a 
singularity and, after expanding (inflating, if u > 0), shrink to another, or starting 
in a non-singular state, collapse to a singularity. I also conclude that the general 
Bianchi- VII;, (with h 7^ 0) models show strong curvature singularities producing a 
complete collapse of the homogeinity surfaces to 2D-manifolds, to ID-manifolds or 
to single points. Our analysis depends crucially on the introduction of the so-called 
intrinsic time, as the product of the JBD scalar field <f> times a mean scale factor 
a 3 = aia2d3, in which the finite-cosmological-time evolution of this universe unfolds 
into an infinite $-range. These universes isotropize from an anisotropic initial state, 
thence I conclude that they are stable against anisotropic perturbations. 
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1. Introduction 

In the Jordan-Brans-Dicke scalar-tensor theory of gravity (JBD) (Jor- 
dan 1959, Brans and Dicke 1961) a massless scalar </>-field is introduced 
in addition to the pseudo-Riemannian metric of the spacetime V4 oc- 
curring in Einstein general relativity (GR); supposedly, this long range 
field is generated by the whole of matter in the Universe according 
to Mach's principle (Dicke 1964); furthermore, the inclusion of this 
field allows Dirac's idea about the secular variation of the gravitational 
constant G to occur in JBD cosmology through <p ~ 1/G. The scalar 
field in JBD acts as an additional effective source of V4 geometry, and 
it is coupled to the tensorial degrees of freedom of the theory by a 
constant parameter co (Ruban and Finkelstein 1975). The value of to 
can be estimated from astronomical observations as |cj| ~ 500 to be 
in accord with current observations. The theory, however, does not 
seem to predict anything different from GR with enough supporting 
observational evidence. 
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However and in spite of what I have just said, there is a great deal 
of interest in JBD (and other scalar tensor theories) due to extended 
inflation and pre-big bang ideas where scalar fields can solve some of 
the problems of inflation (Steinhard 1993, Gasperini and Veneziano 
1994) and may, even, be considered the cause of it (Cervantes-Cota 
and Chauvet 1999). Such scalar fiels can be easily identified with the 
JBD </>-field; moreover, since the inception of superstring theories which 
lead naturally to a dilaton theory of gravity with mandatory scalar 
fields, the importance of JBD theory has been greatly increased. This 
comes about since the JBD action functional already includes a string 
sector where the dilaton field (/>£> can be suitably related with the JBD 
scalar field as (f> oc exp(— 0£>). The important role of the 0-field of JBD 
would especially occur at the strongly relativistic stages of the Universe 
expansion in view of the important role that the so-called scale factor 
duality plays in the pre-big bang scenario (Gasperini and Veneziano 
1994, Clancy et al. 1998). Therefore the importance of studying the 
evolution of the homogeneous and anisotropic Bianchi universes in JBD 
cosmology, since the early cosmological expansion of the Universe can 
be determined by an anisotropic but homogeneous vacuum stage with 
a non vanishing scalar field. The importance of the scalar field can be 
ascertained, for example, by the isotropization it is known to cause in 
the homogeneous but supposedly anisotropic Bianchi universes (Nuhez- 
Yepez 1999, Chauvet and Cervantes-Cota 1995). On the other hand, 
results obtained in JBD models can be interpreted as induced gravity 
cosmological results corresponding to very early epochs of the universe 
by identifying the JBD scalar field with 2tt&$>/uj where 3> is the SU(5) 
isotensorial Higgs field and to is the JBD coupling constant and recur- 
ring to a spontaneous symmetry breaking mechanism (Cervantes-Cota 
and Chauvet 1999). 

In this work I analyse the early epochs of a supposedly anisotropic 
Bianchi type VII vacuum universes and show that these exhibit curva- 
ture singularities which make their surfaces of homogeinity collapse to 
2D-manifolds, to string-like ID-manifolds or, even, to single points. Fur- 
thermore, I show that despite the suppposed anisotropic behaviour of 
the Bianchi type universes, a Bianchi VIIo universe evolves isotropically 
and behaves from the start as a Friedman-Roberson- Walker (FRW) 
universe in spite of possible anisotropic perturbations. The initial an 
final states of this VIIo universe depend on the value of the coupling 
parameter, uj, of the JBD theory; the evolution starts from a singu- 
larity, expands and then reaches another singularity, for to > 0, the 
expanding phase is inflationary. The coordinates used to perform the 
analysis, which includes the scaled (by the local volume on the surfaces 
of homogeneity) 0-field as an intrinsic time coordinate, <!>, are useful for 
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studying the model evolution until the singularity is reached (Nuhez- 
Yepez 1999). Furthermore, contrarywise to the predictions of General 
Relativity (GR), this class of cosmological models isotropize as a result 
of the interaction with the JBD-scalar field not mattering the initial 
values of the scale factors, therefore, I can say that they are stable 
under general anisotropic perturbations. 



2. Bianchi-type VII field equations 

Let us write for the line element of the spacetime, using signature +2 
and natural units c = G = 1, in the so-called sychronous coordinates 
as 

ds 2 = -dt 2 + hijit)^^, (1) 

where the hij(t) is the metric on the surface of homogeinity assumed 
to depend only on t, the synchronous or cosmological time, u l are the 
one-forms (Ryan and Shepley 1975) expressing the properties of the 
3-sur faces of homogeneity whose specific values, appropriate for the 
homogeneous but anisotropic Bianchi VII model, are 

cj 1 = ai ([t] — ku]dy — vdz) , (2a) 

cj 2 = ci2 {ydy — [r] + kv\dz) , (2b) 

cj 3 = a^dx, (2c) 

u A = dt, (2d) 

where rj = exp(— kx) cos(Mx), v = (— M _1 ) exp(— kx) sin(Ma;), k = 
h/2 and M = (1 — k 2 ) 1 ! 2 . Inserting the line element (1), with the forms 
(2), into the JBD vacuum field equations (Ruban and Finkelstein 1975), 
I get 

- 2 (ln ai ) + -(lnaO-Ona^as) + ^(Ina.) ^- j + 

+ Aia{ 2 + Sifo + Fife = 0, i = 1,2,3. (3) 
I additionally have the so-called constriction equation, 



Jt {lnai) Jt {lUa2) + l (lnai) l (lna3) + Jt {lUa2) Jt 0na3) 



d . / 6 \ lo 



2 



+ — (In 010203) I T I - 2 I i J + A ^ + E ^ 2 + = °'( 4) 
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this is an equation of the Raychaudhuri type (Raychaudhuri and Modak 
1988, Chauvet et al. 1992); finally, the scalar field comply with 

4 ( [010203]^ ) = 0, (5) 



dt\ L 1 dt r 

where = (ai/(2a,ja,k)) 2 and the indexes k are to be taken in cyclic 
order of 1,2,3. Equations (3), (4) and (5) are written in the standard 
form we introduced for solving the Bianchi models in JBD (Chauvet 
et al. 1992); the values for the constants appearing in them are: Ai = 
4M 2 - (5/2), A 2 =Ai- 2, ^ 3 = 0,Aa = Ai- 1, £i= £3 = ^=^2 = 
-2, £ 2 = T 3 = 2, £ 4 = F 4 = -1, and M 2 = 1 - {h 2 /A). Notice that 
according to these relationships, h must be restricted to be \h\ < 2 
(this can also be seen from the differential forms in (2)). Also, the 
equations for the Bianchi VIIo model can be particularized from the 
general ones, just by taking h = 0. Given certain features of the metric 
and the concourse of the field equations, the model leds to the following 
two additional relationships 

|(ln ai )-|(lna 2 )=0, (6) 



and 



, . : 0. (7) 

Equation (6) always provides an additional relationship between the 
two scale factors, a\ and 02, not mattering what the value of the 
/i-parameter; on the other hand, equation (7) becomes just a trivial 
identity when h = not restricting in any way the values of the scale 
factors. The case h 7^ and some of its consequences are analysed in 
subsection 3.2. 

2.1. Scaling the </>-field 
From equation (5), I easily obtain 

aia 2 a 3 — = <p , (8) 
at 

where <pQ is an integration constant; thus I can introduce the scaled 
0-field, <!>, as <& = (f)/(f>o, which coincides with the so-called intrinsic 
time (Nunez- Yepez 1995). From equation (8), I easily get 

dt = (a 1 a 2 a 3 y 1 d$. (9) 
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This equations shows that $ is a monotonic function of the synchronous 
time t, $ can hence be used also as a time reparametrization useful for 
solving equations (3). In fact, $ has been found useful for analysing 
the Bianchi vacuum models in several situations (Chauvet et al. 1991, 
1992, Carretero-Gonzalez et al. 1994, Nunez- Yepez 1995). For the sake 
of convenience, let us introduce the notation ( )' = <9$ and, defining 
the Hubble expansion rates as Hi = (lnaj)', the reparametrized field 
equations become 

H[+ -£+ Ji<4 + Ki4 + Malat = 0, i = 1,2,3 (10) 
and the constriction equation becomes 

$ 2 <I> 2 

+ J 4 a\ + /C 4 a| + A/4 = 0, (11) 

the specific values for the constants appearing in equations (10) and 
(11) are combinations of the constants previously used: Ji = Jz = 
£1 = /C 2 = - 1/2, J 2 = £3 = 1/2, Ji = /C 4 = - 1/4, M = 4M 2 - (5/2) , 
A/2 = A/i — 2, A3 = 0, A/4 = Ai — 1. The specific form chosen to write 
the equations and the parameters just emphasizes the relationship with 
our previous work (Chauvet et al. 1991, 1992, Nunez- Yepez 1995, 1999) 
on exact solutions for vacuum Bianchi models in JBD. 



3. Solutions for the vacuum Bianchi type VII universes 

For solving the equations of the Bianchi-type VII model, I found con- 
venient to address separately the specific VIIo and the generic VII/ t 
Bianchi models. The exact solutions of the next subsections are ob- 
tained using essentially the method used in (Chauvet et al. 1992). 

3.1. Bianchi type VII 

In the reparametrized formulation of the equations for the anisotropic 
homogeneous metric of Bianchi type VIIo, solutions can be obtained 
for the case of a Bianchi-type VIIo; the specific solution depends on the 
sign of the quantity 

A = -4(8 + 1/4), (12) 

where 
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B = &HI - -H X - —H[, (13) 

is a constant, i.e. is a first integral of the scaled system (3) and (4), 
that depends on the Hubble expansion rates. In this way I find that out 
of the possible solutions of equations (10), the only physically plausible 
is the one corresponding to the case A < (the A > or A = cases 
can be shown to led to negative or even complex scale factors). For 
some details see (Nunez- Yepez 1995). The only physically admissible 
solution can be explicitly written as 



«i($) = (^^) 1/4 ($cosh [-V0B + 1/4) hi(/<I> 2 )]) 



-1/2 



, (14) 



where cq and / are positive integration constants. The other two scale 
factors can be easily obtained from ai(3>), as follows from (6) and (10), 
they are 



a 2 ($) = c ai($), (15) 

a 3 ($)= 2- 1 / 2 Co a 1 ($). (16) 

These results explicitly show that the three scale factors are propor- 
tional to each other; this means that the Bianchi-VTIo model, despite 
what I could have anticipated, shows an isotropic evolution; it also 
implies that the shear, vorticity and acceleration of the reference con- 
gruence all vanish. The vacuum Bianchi-VTIo JBD universe thus be- 
haves as a Friedmann-Roberson- Walker (FRW) space-time — in a way, 
this is not totally surprising since Bianchi cosmologies correspond to 
the simplest deviations from a FRW environment. Figure 1 shows the 
evolution of the scale factors of the model as a function of the intrinsic 
time (or scaled 0-field) <J>, for the values u> = 48,0,-4/3. The local 
volume on the surface of homogeneity is then 

V = ^f. (17) 

In the case go > the model exhibits inflationary behaviour, as 
can be easily proved from equations (14)— (16), reaches a maximum 
volume in a finite cosmological time and, in the process, erases any hint 
of the primordial anisotropies. It then shrinks again to a singularity 
(Nunez- Yepez 1999). The constriction equation implies the following 
relationship in our case 
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(18) 



which directly implies that B = ui/ '6; from here I notice that to have 
meaningful solutions the coupling parameter has to be restricted to 
uj > —3/2 or, in terms of B, the restriction is B > —1/4. The con- 
striction equation is found thus to depend only on uj, whose specific 
value is enough to determine the evolution of the Hubble expansion 
rates through the nonlinear but elementary equation (13) — very similar 
equations are valid in all the integrable Bianchi models (Nunez- Yepez 
1995). I have also integrated the equation obtaining (in the case to > 0) 



Hi($) = v^Btanh ( arctanh(if $o/v / ^) + 2v / £log($ /$) J /$; (19) 



notice that, as the scale factors are proportional to each other, the 
Hubble rates are all the same at all times: H\ = H2 = H%. From the 
solution it can be seen that the Hubble expansion rates, after reaching 
a minimum value, vanish asymptotically as ~ log (<£)/<!> when <3? — ► 00. 

Using equations (9) and (17) I can obtain the dependence of <!> on 
t, as follows 



This integral can be explicitly evaluated in terms of the Barnes ex- 
tended hypergeometric function (Rainville 1980) but, in practical terms, 
is better to calculate it numerically. It should be noticed that on inte- 
grating (20) from to 00, the cosmological time reaches a ^-dependent 
finite value: t e (Nunez- Yepez 1999); I cannot reach beyond this specific 
value of the cosmological time in our approach. Though I can obtain t 
as a function of <I>, I cannot invert it to obtain explicitly $ as a function 
of t. Nevertheless, figure 2 exhibits the dependence of the field <I> on 
t, showing the enormous change that occurs in over a very small 
span of t values. This also shows that, asymptotically, $ grows without 
bound whereas t approaches a certain finite value i e (— 3.624 in this 
specific case); on the other hand as t — > 0, <I> vanishes. Notice that 
these conclusions can be applied with no changes to the scalar field 
4>, excepting when 4>q = 0. This behaviour has been instrumental for 
explaining, in a Bianchi IX cosmology, the vanishing of the maximum 
Lyapunov exponent calculated in the intrinsic time $ in spite of it being 
positive in the synchronous time t (Carretero-Gonzalez et al. 1994). 

I can now relate the behaviour of $ with that of 01; since, as in figure 
1 shows, the universe can begin with a singularity and then, as <3? begins 




(20) 
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grows, the universe expands as ~ <3? a , where a = + 1/4 — 3/2. But 
this expansion from a singularity only occurs if B > (hence, only if 
<jj < 0), otherwise the universe can start from a non singular state 
and proceed to collapse to a singularity. In the expansion case, the 
universe rapidly reaches a maximum volume F m ax, that can be easily 
calculated from (14) and (17), and then it shrinks until it reaches again 
a singularity; see figures 1 and 3 and the comments in section 4. As a 
conclusion of the previous discussion I can say that our chosen time 
coordinate, is seem to be appropriate for studying the early stages 
of the evolution and certain features of the singularity (see section 4 
and figure 2) however, it does not allow us to analyse the behaviour 
beyond it. However, as I am mainly interested in the early epochs of 
the universe, this is of no importance. 

3.2. BlANCHI TYPE VII/j 

From equations (6) and (7), a few relations can be obtained for the 
scale factors irrespective of the value of h / 0; from (6) I get (15) (the 
same equation than in the case h = 0). From equation (7), I might get 
a2 = 0, implying that a± = and that 03 can take any value; then 
the model collapses into essentially a spatial one-dimensional manifold: 
i.e. a string- like object. Other possible options allowed by equation 
(7), imply that 02 G13 —>■ 00, thus the spatial 3-surfaces of homogeinity 
are seen to collapse into 2-surfaces. The important point here is that 
for any choice of values for the scale factors in equation (7), the model 
is found always to spatially collapse. Every Bianchi-type VII/j (with 
h 7^ 0) universe is thus highly singular (see section 4). 

As I have been able to get the scale factors for the especific cases 
addressed in this paper, I have obtained exact vacuum solutions for 
the Bianchi VII models. The important conclusion is that all solutions 
found for the generic Bianchi VII model with h = 0, and h 7^ — with 
the restriction \h\ < 2 — describe collapsing universes. The behaviour 
of the scale factors in the generic cases (with h = 0) is exhibitted in 
figure 1. 



4. Curvature singularities in the Bianchi-type VII universes 

In this section I study the curvature singularities present on the Bianchi- 
VII universes, though some very specific singularities were discussed in 
subsection 3.2. 

I say that a universe is singular if the value of the Ricci scalar R = 
g ab R a b along a congruence of geodesies diverges, \R\ — > 00, whereas the 
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associated affine parameter tends to a finite value (Wald 1984). If this 
happens, then I say that there exists a curvature singularity and thus 
that I am dealing with a singular universe model. Notice also that, due 
to the choice I made for the uj parameter, here the singularity occurs 
when R — > — oo (figure 3) rather than the other way round. 

The problem is basically how to choose an appropriate geodesic 
congruence. For our Bianchi VII models, I have chosen as the proper 
congruences the world lines of test observers (time-like geodesies) whose 
affine parameter is the synchronous time t. Thence, the scalar curvature 
can readily be shown to be 



Using only equations (5), (10) and (11) I can rewrite R for the vacuum 
Bianchi- VII model in terms of the coupling parameter oj and the scalar 
field 4>, as follows 



it is important to notice that expressions (20) and (21) do not depend 
on the h- value in any way, they are valid for all the Bianchi-type models 
and not only for the type VII I am discussing in this paper (Niinez- 
Yepez 1995). I have expressed the scalar curvature in two different 
ways, equations (20) or (21), both are important because they exhibit 
the explicit dependence of R on the scale factors and <p and its time 
derivatives, or on the coupling paramenter to and the scale factors (21) 
and, besides, they have both a certain simplicity. For R — > oo in (21), 
all that is needed is that at least one of the scale factors (ai, a2, or 
03) or the reescaled scalar field <£, or just the scalar field (ft, vanish at a 
finite value of the synchronous time t. Notice also that I can regard the 
evolution of the Bianchi-type VII universe — equations (3) and (4) — as 
driven by the curvature R, this is especially true near the singularity. 

Notice that the plot of t against in figure 2, which does not depend 
on h in any way, shows the rather small range of <I>-values in which the 
expansion of the model universe VIIo occurs, as I can see comparing 
with figure 1; this also corresponds to the region free of singularities 
in the model (VIIo), as can be seen on comparing with figure 3. On 
comparing figures 1 and 2, the role of the curvature scalar in governing 
the expansion can be qualitatively described as follows: a strong cur- 
vature prevents expansion; it is only when curvature is small that the 
full extent of expansion is reached but, as soon as the curvature is large 




(21) 




(22) 



asslatex.tex; 7/02/2008; 1:47; p. 9 



10 



H. N. Nunez- Yepez 



in magnitude again, contraction sets in. This qualitative behaviour is 
in accord with current ideas (Turner and Weinberg 1997, Clancy et al. 
1998). 

In the general VII/j JBD model, the universe is always singular, as 
can be easily seen from equation (7) and the results of section 3.2; in 
this sense, I say that the model VII^ is completely singular. 

5. Concluding remarks 

The supposed homogeneous and anisotropic Bianchi VII model shows 
in fact an isotropic expansion in the case in which h = 0. On the 
other hand, this article shows that the dynamics of the early stages 
of the expansion in the same JBD model depends on just one of the 
scale factors (that, here, I choose as a±) and that much of the universe 
evolution depends on the sign of uo (or on that of B), since equation 
(13) determines the Hubble expansion rates. This is as happens in 
other Bianchi models (Chauvet et al. 1992). I have also obtained the 
dependence of the three scale factors on the reescaled field <!>. I 
have concluded that the three scale factors are proportional to each 
other thence the expansion is not anisotropic, as it is usually assumed 
to be in this model. In fact, I have shown that a Bianchi- VIIo JBD 
vacuum universe is basically equivalent to a FRW-spacetime with no 
shear, no rotation, and no acceleration; these are characteristics of our 
actual Universe which can be inferred from the smoothness (and other 
observational evidence) of the cosmic microwave background radiation. 
An important point of the analysis is that it shows that, even starting 
with supposedly anisotropic models, the inclusion of a scalar field (the 
JBD-field) can drastically isotropize the behaviour and thus offers the 
possibility of coordinating an arbitrary anisotropic vacuum initial state 
with the observed isotropic properties of the actual Universe — in spite 
of the issue of the primordial global anisotropy (Christoulakis et al. 
2000). This property shows also that the universe dynamics is stable 
against anisotropic perturbations. Moreover, this feature does not seem 
to depend on u in any way, it could be valid even with a dynamic 
coupling parameter as needed for extended inflation models. 

The analysis is made possible by the choice of the intrinsic time <3? 
as a coordinate. It should be clear, however, that from such coordinate 
patch I cannot obtain information about what happens with the model 
in the far future (i.e. what happens as t — > oo); from this point of 
view, the analysis performed refers to the early stages of the model 
evolution. I have shown that any primordial anisotropy can be erased in 
a finite cosmological time due to the effect of a scalar field. On the other 
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hand I have also shown that certain Bianchi models in the JBD theory 
can be completely singular, as the type VII— h (h 7^ 0) analysed here. 
Notice also that the models show inflationary behaviour when u > 
and that any primordial anisotropy is smoothed and eventually erased. 
A possible realistic setting for our model could be that a primordial 
vacuum anisotropic state of the universe that, due to the presence of 
a scalar field, appears isotropic and essentially FRW after, and even 
during, inflation. 

Nevertheless, for the case of the VII/i, h ^ 0, vacuum model I can 
do conclude that, not mattering the choice made in equation (7), the 
universe always collapses. According to the discussion in section 3.2, the 
homogeneity 3-surfaces of the universe spatially collapse into 2-surfaces 
or into one dimensional objects or, even, into a single point; this model 
universe always collapses to a permanent singularity. Such behaviour 
may make this model of certain interest if one is just interested in 
singular behaviour. In this respect, the dependence of the universe 
dynamics on the curvature scalar is worth pinpointing. 
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Figure Captions. 

Figure 1. 

The behaviour of the scale factors in a Bianchi-type VIIo universe 
as a function of the intrinsic time $ is shown. The three scale factors 
are proportional to each other. I have used the values / = 1, cq = 2, 
B = w/6, for the parameters in equations (14)-(16). 

a) The case oj > 0, specifically I took oj = 48. 

b) The case uj = 0. 

c) The case uj < 0, specifically uj = —4/3. 
Figure 2. 

The relationship between the cosmological time t and the intrinsic 
time <3? showing a monotonic behaviour in every case and the finite 
range of t values associated with an infinite range in <E> values. Notice 
the short span of $-values where the main changes in the universe take 
place. 

a) The case uj > 0, specifically uj = 48. 

b) The case uj = 0. 

c) The case uj < 0, specifically uj = —4/3. 
Figure 3. 

The absolute value of the scalar curvature R in the model dealt 
with here is shown in the range < R < oo. I am graphing arctan(|i?|) 
against arctan(<J>). In this way, small t- values correspond to small <J>- 
values, but large values of $ correspond to t ~ t e . The graph was 
compactified in such a way that a in the horizontal axis is really 0, 
but when it/2 appears, it should be interpreted as oo. The three cases 
a), b) and c), shown are the same as in figures 1 and 2. 
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